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Abstract 

The non-equilibrium two loop self-energy is reexamined in the frame- 
work of a scalar quark and gluon model, with specific attention to 
terms which do not give rise to the standard two — * two particle col- 
lisional terms in the semi-classical Boltzmann equation. It is shown 
that most of these terms contribute to renormalization of component 
lines at a lower level, rendering the theory correct to g 4 m 4 . This result 
can be generalized to a higher number of loops. The remaining terms, 
which do not fall into any physical category, are shown explicitly to 
vanish. We then examine the possibility that pinch singularities could 
arise in this theory, and demonstrate that this is not so for the case 
of equilibrium and small deviations from equilibrium, at the two loop 
level. 
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I. INTRODUCTION 



Due to the short time span of a heavy-ion collision, it is believed that a non- 
equilibrium description of such a collision may be necessary to give a complete and 
adequate description of the particle evolution. Several reviews of the non-equilibrium 
approach in general may be found in M. For heavy ion collisions, heuristic models 
have been developed, which are guided by the approach of Boltzmann, that presumes 
that only two body collisions can occur. These can then be generalized in an heuris- 
tic fashion |2| . On the other hand, from the field theoretical point of view, there have 
been many attempts to generate the correct and appropriate non-equilibrium the- 
ory for a given Lagrangian (see, for example || . Often these stop at the mean-field 
level, due to the inherent complexity of the non-equilibrium formalism per se. In 
some previous works, the two-body collisional term has been calculated for a static 
potential interaction, and under many simplifying assumptions in [[|. Although a 
graphical series was given for the two loop self-energy following from a Walecka type 
interaction in || , the physical interpretation of all the terms that occur in this series 
or direct expressions for them was not given. Note that the difficulty lies in the fact 
that the non-equilibrium formalism generates far more diagrams of a given topol- 
ogy than does the equilibrium zero temperature or Matsubara finite temperature 
formalism which generates Feynman graphs in the standard way. In a broad sense, 
the ultimate goal which we have, is to obtain a complete physical understanding of 
all the additional terms (over and above those which are expected from a Feynman 
diagram point of view) that are generated in the non-equilibrium formalism, in order 
that a precise formulation in, for example, the study of non-equilibrium evolution 
of distribution functions in QCD, rather than simply a heuristic approach, can be 
made. At this stage, however, we restrict ourselves to far simpler systems. 

In this paper, we examine the collisional integral for a scalar quark-gluon model 
in the Hartree-Fock approximation and for the two loop diagrams in greater de- 
tail than before |J, placing specific emphasis on the extra terms that are naturally 
generated in the non-equilibrium formalism. In this study, (a) we are able to show 
that two-loop diagrams that were previously ignored as they do not contribute 
directly to the two-particle scattering terms expected from classical Boltzmann the- 
ory, do in fact have an important physical relevance: all of these remaining graphs 



2 



are required to correct the lower order tree level graphs which give rise to absorption 
or radiation processes and which also occur as matrix elements in the Boltzmann 
equation to one loop, implying that a simulation of a non-equilibrium process with 
this or similar models requires inclusion of such terms, (b) Our study is carried out 
within the Schwinger-Keldysh formalism J7[ for non-equilibrium processes. Within 
this (and other) formulations of real-time finite temperature processes ||, the issue 
of pinch singularities arises (see for example P-pH). This addresses the question 
as to whether the theory is well-defined or not, in that it possesses divergences 
that arise apparently from integration along the real axis of products of retarded 
and advanced Green functions GrGa- Such products occur naturally in the so- 
called 'retarded-advanced' (RA) version of the non-equilibrium formalism, which is 
an alternative way of presenting the Schwinger-Keldysh Green functions [0. Alter- 
nately, pinch singularities manifest themselves as products of delta functions within 
the standard representation which we follow here. Within our model, we investigate 
the occurrence of pinch singularities in the series for the self-energy up to two loops 
and come to the conclusion that no pinch singularities occur. 



To be specific, we utilize the scalar partonic model fI5,n|, defined by the La- 
grangian 

-gm^\T a ){(T r )r<P hmXa ,r - ^UJrstX^X^ (1-1) 

in which scalar quark fields <fi 1 ' 1 are coupled to scalar gluonic fields x a,r Yl& a cubic 
interaction. The color indices i,l = l..N c , while a,r = l..(N c — 1), reflects the 
product groups SU(N C ) x SU(N C ). 



II. COLLISION TERM TO TWO LOOPS 

In order to recover the semi-classical extension of the Boltzmann equation, it is 
necessary to evaluate constructs of the form |||| 

J coll = Tl~ + (X,p)D + -(X lP ) - U + -(X,p)D-+(X,p) 

_ rgain _ rloss (O 1\ 

— J coll J coll ) y^- 1 ) 
where D^(X,p), i = +, — are Green functions of the standard Schwinger-Keldysh 
kind, and which are given in Appendix A. The transport and constraint equations 
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are briefly summarized in this Appendix also. The IP- 7 are the Wigner transformed 
self-energies of the associated particle. To be specific, we consider for the moment 
the self-energies and Green functions pertaining to the quark sector, U^(X,p) = 
E u (X, p) and D l i (X, p) = S 1 ^ (X, p) . The analysis for the gluonic sector follows then 
simply by inspection. Let us consider first the loss term. This means we need to 
calculate E + ~(X, p). The quark self-energy to one loop is shown in Fig. la) which is 
the Fock term, and to two loops in Fig.2. As there are many graphs at the two loop 
level, we make a structural classification of these according to their topology. We 
identify rainbow graphs (R), ladder graphs (L), cloud graphs (C), exchange graphs 
(E) and quark-loop graphs (QL). In |J, we were able to show explicitly that the 
Fock term of Fig. la) leads to the requirement that the process qq — > g, as given 
by Fig. lb), be incorporated into the associated Boltzmann equation at this level. 
We have also shown that the seven diagrams R a), La), C a), C b), E a), E b) 
and QL a) are all essential for constructing the direct cross-sections of the processes 
qq — ► qq, qq — > qq, qq — > gg, and qg — > qg which are shown in Fig3, and which then 
appear in the newly calculated Boltzmann equation that should be correct to the 
two loop level. Two questions now arise naturally. Firstly, one may ask what the 
purpose of the remaining thirteen graphs is, and secondly, more subtly, would they 
harbour divergences such as arise through pinch singularities, thereby rendering the 
theory ill-defined. In a previous paper, the remaining thirteen two-loop self-energy 
diagrams shown in Fig.2 were argued heuristically to be 'non-leading', and graphs 
of this type have up to now simply been ignored. Here, we wish to demonstrate 
precisely that these graphs in fact either vanish or contribute to corrections of order 
g 3 m 3 to the lower order process qq — > g. To demonstrate this, we arbitrarily examine 
the set of quark-loop diagrams. The QL a) graph of Fig.2 leads directly to the qq 
and qq cross sections of Fig. 3. On the other hand, the three graphs, QL b)-d) of 
Fig.2 which were not required for the evaluation of the cross-sections as derived for 
Fig.3. Let us examine each of these terms in turn. We commence with the diagram 
QL b) of Fig.2 which is given by 

V,(2)+- (v \ 4 4772 f d4 Pl d% d% d% s4r(4)/ \ 

£ q -Loop, b) (^) = -gm F ^ op Jj^j^j^j^W * } (p-pi - ft) 

x(2tt) 4 <^(p 2 -p 3 +p 4 ) S+-{X, Pl ) G + -(X,p 2 ) 

xG ++ (X,p 2 ) S ++ (X,p 3 ) S ++ (X,p,), (2.2) 
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where -Fq.Loop i s the color factor for the two color SU (N c ) groups. The corresponding 
loss term of the collision integral to this self-energy is defined by 

^S° q -Loo P ,b) = -vCoop,^ 1 -/ = E p ,p)f g (X,p). (2.3) 

In this expression, the product iS + ~(X, pi)iG + ~(X, p2)f q (X, p) occurs. Inserting 
the quasiparticle approximation for the Green functions that are given in appendix 
A, we obtain for this product the sum of four terms: 



T x = 


7T 

E~i 


^ r 5(E 1 -p 1 )5(E 2 -p° 2 )f q (X,p 1 )f 9 (X,p 2 )f q (X,p) 




T 2 = 


7T 

E~i 


^-6(E 1 -p 1 )S(E 2 +p° 2 )f q (X,p 1 )f g (X, -p 2 )f q (X,ft) 

Hi 2 




T 3 = 


7T 

E~i 


^-S(E 1 +p° l )S(E 2 -p 2 )f q (X, -pi)f g (X,p 2 )f q (X,$ 

Hi 2 




T 4 = 


7T 

E~i 


^ r S(E 1 +p 1 )5(E 2 +p° 2 )f g (X, - Pl )f g (X, -p 2 )f g (X,p). 

M, 2 


(2.4) 



By attributing / to incoming particles and / to outgoing ones, we see that Ti...T 4 
correspond to the processes q — > qg, qg — > q, qq — > g and qqg — >0. Since the 
quarks are massless while the gluons are endowed with a finite mass, the processes 
corresponding to T 1; T 2 and T 4 are kinematically forbidden as also occurred in the 
discussion of the Fock term in ||. One thus has one remaining non- vanishing con- 
tribution iS + ~(X, pi)iG + ~(X, p 2 )f q (X, p) = T 3 . This product is now inserted into 
Eq. ( |2.3D . In the resulting expression, we interchange p\ with —p\ and on performing 
the Pi, p 2 and p 4 integrations, we find the result 

J coll,q-Loop,b) " * Ep 9m Vloo P J (3^33^ ^Tt^E, (2tt)* ( } ^ + Pi W 

xG ++ (X,p 2 )5 ++ (X,p3)5 ++ (X,p3-p 2 )/ 9 -(X,^)/,(X,p 2 )/ g (X,p) 

(2.5) 

as the remaining contribution of the QL b) graph to the collision integral. 

Now, in Fig.4, we independently examine the Feynman diagrams for the process 
qq — > g up to order g 3 m 3 . The scattering amplitude associated with Fig. 4a) has 
purely a point-like structure with color groups occurring: 

-%M%^ g = -igmt%®t\ m (2.6) 

while from Fig. 4b), one has 
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- iArfJU = 9 3 rn 3 [t*Mt b t a )} ® [t^tr(tT)]G— 

where is the matrix of the color group in the representation of the quarks. Now 
note that using the fact that [iG ]> = iG ++ and -F q _ Lo op = t^t b j i ti{t b t a ), one can 
rewrite Eq.( |2.5| ) in terms of these matrix elements, i.e. 

E p J (2n) 3 2E 1 (27r) 3 2E 2 ' 

xM^glM^f^X^f^X^fgiX,^), {21 



t(2)1oss 71 / UjPI «7>2 /„ Ur(4)/ , X 



illustrating that the cross term between these two processes, denoted symbolically 
as a)b, is derived from the self-energy diagram QL b) of Fig. 2. The gain term can 
be obtained by replacing / with / and vice versa in Eq. fl2.8| ). 

In a similar fashion, the collision integral can be constructed from the quark- 
loop diagram QL c) in Fig.2. One obtains an expression for the loss term as in 
Eq.( [2.5|) with iG ++ iS ++ iS ++ replaced by the combination iG iS iS . Again 
^coii'q-Loop b) can t> e expressed by the scattering amplitudes of Eq. ( |2.6| ) and (|2.7|) : 



7 (2)loss VT f (Ppi <fip 2 (n Ur(4)/ , \ 

x[M^M%^ 9 f q {X^fq{X,p 1 )f g {X,p 2 ) : (2.9) 

i.e. the second cross term tfa required in building a cross-section of the basic com- 
ponent a) and b) of Fig. 4 is required. 

In an analogous fashion, one can show that the rainbow diagrams R b) and c) 
lead to a collision integral containing Aifq^ tg [A / lqq-*g\ ji and the hermitian conjugate 
of this product, the ladder diagrams Lb) and c) to a collision integral containing 
Aiql^glAi^q^g}^ and its hermitian conjugate, the cloud diagrams C c) and d) to a 
collision integral containing A4ql^ g [A4ql^ g }^ and its hermitian conjugate, and finally 
the exchange diagrams E c) and d) to a collision integral containing M. a q \^ g {M e q \^ g \ 1 '' 
and its hermitian conjugate. Note that in this fashion, we are able to account for all 
mixed diagrams that would occur in the construction of the \M. q q^ g \ 2 up to order 
g 4 m 4 , with the exception of the diagram g) of Fig.4. This graph does not enter into 
the collision integral, as it is a renormalization diagram for the incoming quark, for 
which the momentum p is fixed externally. 
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Returning to our explicit example of the quark-loop self-energy of Fig. 2, one 
sees that a simple graphical interpretation can be applied to each figure. A rule in 
which all lines that are connected by ± and =F are cut in a single path, separates the 
graphs QL a)-c) into their component matrix elements. This is illustrated in Fig.5. 
This procedure, however, cannot be applied uniquely to the graph QL d), nor for 
that matter to the remaining graphs which are not required for construction of the 
mixed terms or direct contributions to the cross-sections, i.e. the graphs R d) and 
L d). We are thus now left with the three graphs QL d), R d) and L d) which at 
first sight fit into no apparent scheme, and which therefore may present difficulties. 

We commence with the investigation of QL d). For this self-energy we obtain an 
expression as in Eq.( [2.2|) with the product of the five Green functions replaced by 
S + ~ (X , pi)[G + ~ (X , p 2 )} 2 S~ + (X , p 3 ) S + ~ (X ', p 4 ) . In the quasi-particle approximation 
the off-diagonal Green functions are on mass shell: 

P 2 i=pi=pl = (2-10) 
p 2 2 = m 2 (2.11) 

In addition to this, the two 5-functions for the energy-momentum conservation of 



Eq^^) have to be fullfilled. Therefore we can write for example 

= vl = (pa + Pa? = 2P2P4 + m 2 . (2.12) 

Together with Eq. ( [2.1 1| ) this leads to 

1 11 
P4 = --P2 => p\ = = I™ 2 > 0. (2.13) 

This is inconsistent with Eq.( [2.10 ) and therefore ^^XoopA) nas ^° van i sn - A similar 
investigation of the graphs R d) and L d) shows that they vanish for the same reason. 

So far we have derived the results for two loop self-energy diagrams in detail. 
Let us summarize now the results for self-energy diagrams with more than two 
loops. The result obtained in this section can be generalized easily to three and 
more loops. Including the three loop self-energy into the collision term leads to 
cross sections of all possible processes with two (three) partons in the initial state 
and three (two) partons in the final state. Furthermore, they give all additional 
corrections of order g 6 m 6 to the process qq — > g, e.g. the product of scattering 
amplitudes of the diagrams b)-f) in Fig. 4, and also corrections to all scattering 
processes of two partons into two partons, qq — > qq, qq — > qq, qq — > gg and qg — >■ qg. 
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III. PINCH SINGULARITIES 



To investigate the issue of pinch singularities, let us make two simplifmg as- 
sumptions: a) the self-energies are calculated at zero momentum and b) we do not 
distinguish between quarks and gluons. The latter assumption reduces the number 
of generic two loop self-energy diagrams from five to two, shown in Fig. 6. The first 
type, the rainbow diagram, is shown for and in more detail in Fig. 7. From 
these two self-energies, the retarded self-energy can be constructed as 

E r R = ET + ^R + - (3.1) 

Since the retarded self-energy is a physically relevant property, that enters, for ex- 
ample, into the constraint equation Eq. (|^5|) , we want to investigate the diagrams 
of Fig. 7 in more detail. In each of the diagrams b), c) and f) of Fig. 7, it is possi- 
ble to identify an internal vertex to which lines three off-diagonal Green functions 
are attached. Since these Green functions are on mass shell, this corresponds to a 
decay of an on-shell particle into two on-shell particles of the same species. This is 
a forbidden process, and therefore these three diagrams vanish. Closer inspection 
of Fig. 7 g) shows that the momentum structure of the Green functions is similar to 
that of b), c) and f), and this graph is therefore also vanishing. Consequently, 
is the sum of the remaining diagrams a), d), e) and h): 

d 4 k dH fr „ 



E R (X,0) = (-igm) 4 J 



(2tt) 4 (2tt) 4 l L 
i n — aun-+/ 



-(*- 


-l)D~ 


-(I) 


\k- 


l)D + - 


f (0 


-(*- 


l)D~ 


"(0 


f (A;- 


-l)D + 


+ (0} 



(3.2) 

This construction appears to contain pinch singularities, evidenced by the fact that 
products of D{k) occur, and it is imperative to show that the apparent divergence 
vanishes through cancellation. For simplicity, we consider this problem first in equi- 
librium. We do this for two reasons. Firstly, a correct calculation in equilibrium 
should be entirely free of pinch singularities, as has been demonstrated generally in 



T5|j. In addition, calculations can be simplified dramatically, in such a fashion as 



to allow for an extension to systems slightly off equilibrium in a simple fashion in 
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this direct example. We use the following reasoning: in equilibrium, the Schwinger- 
Keldysh contour is simply one of a family of choices of contour in which the upper 
and lower integration routes are separated by a distance a, and with which a form 
of real time field theory [RTFT] is generated |TH] . In the Schwinger-Keldysh case, 
one makes the choice a = 0. Physical quantites, when calculated, should however 
be independent of the value of a. As it turns out, the particular value a = 1/2, 
which generates the so-called thermo field theory (TFT) formalism |fL7| , gives dis- 
tinct calculational advantages, as products of off-diagonal elements can be seen to 
simply vanish. Given this equivalence, we are free to choose to do our analysis in 
the TFT framework. 

The scalar Green functions, in the TFT formalism, are given as 

(ID— iD~+\ = (Dl + Dp D'p \ 
\iD+- iD++) \ D'p Dl + Dp) 

2nS(k 2 -m 2 ) 9 r/, 2 _ *\ e^oW 2 



/ k 2 -m 2 +ie + w e m \-i 2n5(k 2 m ; e(3 | feo |_ 1 

\?7rA(P e/? ' fc °' /2 -» , 2^{k 2 -m 2 ) 



(3-3) 

■ e /3 1 fc o I — 1 k 2 ~m 2 -ie n e /3 1 feo I — 1 ^ 

so defining Dq, Dq, and D^ 1 . Let us compare these Green functions with the 
ones in Eqs. flAlOD - QAT3I) . In equilibrium, / a (X,p) = f a (X,-p) = l/(e /3 l p °l - 1). If 
we do not distinguish between particles and antiparticles, then the diagonal Green 
functions in the TFT formalism are equal to the ones in Eqs.( |A12| ) and ( |A13| ). The 
off-diagonal Green functions differ from Eqs.( |AlO| ) and ( |A11| ). However, in the self- 
energies considered only the product D~ + D^ occurs, which is the same in both 
types of formalism. This in itself demonstrates that one may use the Green functions 
of the TFT formalism with impunity in the following discussion. We consider first 
the last term of the retarded self-energy in Eq. fl3.2p : 

/ ~B) iD++{k _ i)d++{i) = I { n2 ° {k _ i)d1{i) + Di{k _ i)D(s{i) 

+D p (k - l)D 2 (l) + D p (k - l)D p (l)} 
= J {-Dl(k - l)Dl(l) -Dl(k- l)D p {l) 



1 Our notation introduces a factor i in the first matrix, in order to establish consistency 
with our Schwinger-Keldysh Green functions of Appendix A. This differs from the standard 



use of TFT users [18]. 
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-D fi (k - Z)Dj(f) + Dp(k - ()£>„(<)} 

In the first step, we have made use of the relations ]18 

/ (Sv Dl ° {k - l)Dl ° {l) = ~J ^Y Dl{k ~ l)D * {l) (3 ' 5) 



and 

/ ^ Dl ° {k ~ l)Dp{l) = ~J (^ D2 ° {k ~ l)Dp{l) - (3 - 6) 

In the second step, we have made use of the fact that the fourth term, Dp(k—l)Dp(l), 
has to vanish, since together with D'^k) of Eq . ( fO]) it corresponds again to the decay 
of an on-shell particle into two on-shell particles. With Eq. (|3.4j) , we can rewrite 
Eq.( |3.2| ) for the retarded self-energy as 

E r R (X, 0) = {gmf J ^^D~(k - l)D~(l) 

x {[D—(k)} 3 - 2D—(k)[D-+(k)} 2 + D ++ (k)[D- + (k)} 2 } . (3.7) 
It is now convenient to make use of the following representation of 5-function, 

2n5(k 2 - m 2 ) = — — , - — — 

Ar — m + ze k — m z — ze 

= Dl + Dl (3.8) 

Then the term in curly brackets of Eq. (|3.7|) is evaluated to give 

{...} = [Dl + {Dl + Dl)f B f + [D 2 + ip\ + D 2 )f B ](Dl + Dl) 2 g 2 B 

-2[Dl + {Dl + Dl)f B ]{Dl + D 2 ) 2 g 2 B , (3.9) 

where g B = e /3 l fc °l //2 /(e /3 l fc °l — 1). Using the fact that f B — g B = — f B , we find 

{...} = (Dir + [(Dlf + (D 2 ) 3 ]f B . (3.10) 

This expression is well defined, since no products of D^D 2 occur any longer. With 
the relation 



2 \ dm 2 



d ^ 2 



2 \dm 2 
10 



2n5(k 2 -m 2 ), (3.11) 



we can write the retarded self-energy as 

d 4 k dH 



(2vr) 4 (2vr 



-D-{k-l)D-{l) 



x Ud 



1\3 



1 



d 



2 \ dm 2 



2n5{k 2 -m 2 )f\ 



B 



(3.12) 



We conclude that for the retarded rainbow self-energy no pinch singularities 
occur. 

We now have to consider the cloud self-energy of Fig.6(b). The retarded self- 
energy is constructed in a similar way as for the rainbow diagram. The calculation of 
such a diagram within the framework of 3 theory has been performed in [ IS] . This 
result can be simply taken over for our purposes, and we quote the final expression 
here: 



r d A k d A l 
C K ' ; yy ' J (2tt) 4 (2tt) 4 

x ~D l (k + D p (k - I)] (D l (l)Y + X {1) 



Dhik)) + X (k) 

2 



(3.13) 



where the function 



X (k) = 2DlD p 



n 2 



2P 



k 2 — m 2 



2n5(k 2 - m 2 )f B 



d 



dm 2 



2n5(k 2 - m 2 )f. 



B 



(3.14) 



is free of singularities. Thus we conclude that in equilibrium, the retarded self- 
energies with two loops have no pinch singularities. We conjecture that the inclusion 
of different kinds of particles interacting with each, as in Section [TT], will not display 
pinch singularities - the actual demonstration of this fact, however, is cumbersome. 

Let us now investigate the effect of small deviations from equilibrium, i.e. we 
replace f B = l/(e /3 l fc °l - 1) by f B + 5f B in Eq.Q. In deriving our result for the 
retarded self-energies and Eqs. fl3.12D and (|3.13| ), a central feature is that the 
relationship f B — g 2 B = —f B was used. For small deviations from equilibrium, this 
becomes 



fl 



9 2 b 



/b(1 + /b) 



Ub + SJb] 2 

[f B + Sf B }[l + f B + 5f B } 
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=> /£-</!-> Ub + Sf B ] [f B + 6f B -l-f B - 5f B ] 

= -[f B + Sf B ], (3.15) 

i.e. it is still valid. Therefore we obtain for and Y7 C expressions such as in 
Eqs. (|3.12j) and ( |3.13| ), respectively, with f B replaced by f B +5f B . Consequently, even 



for small deviations from equilibrium, no pinch singularities occur for the retarded 
self-energies up to two loops. 



IV. SUMMARY AND CONCLUSIONS 



In this paper, we have demonstrated that for the two loop self-energy, the graphs 
which are generated in the non-equilibrium Schwinger-Keldysh formalism which do 
not contribute to two — ► two processes are either vanishing or are necessary to renor- 
malize the graphs of one loop order to the same level of the coupling strength g 4 m 4 . 
This applies to all external legs, save the incoming leg, which is fixed externally. 
This result can be generalized to n loops: self-energy graphs containing n-loops lead 
to contributions to the semiclassical Boltzmann equation that are n — > n in nature 
plus include all permutations of the particles on the left and right hand sides 0. In 
addition, however, graphs of lower order are renormalized successively to the same 
order of the coupling constant. We have illustrated that, to the two loop level, 
pinch singularities da not occur in the equilibrium formulation, as expected from 
Ref. [[nj and by a simple extension of the technique here, for small deviations from 
equilibrium. 



V. ACKNOWLEDGMENTS 



One of us (D.S.I.) thanks C. Greiner and S. Leupold for useful comments made 
during the course of this study. 



APPENDIX A: GREEN FUNCTIONS 



In this Appendix, we list the Green functions used in Section [0]. In the 
Schwinger-Keldysh formalism |7j they are defined for the (scalar) quarks (D = S) 

as 
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iS c (x,y) = (T<jt> l (x)<l>V> m (y)) - (0 M (x)) (^ j ' m (y)) = iS~~(x, y) 

iS a (x,y) = (f^ l {x)(j)^ m {y)) - (^{x)) (<j) jj > m (y)) = iS ++ (x,y) 

iS > (x, y) = {f\x)(j)^ m {y)) - (f' l {x)) (f j > m {y)) = iS + ~(x, y) 

iS < (x,y) = ^ m {y)^\x)) - (f>\x)) (<j) jj > m (y)) = iS- + (x,y) (Al) 

and for the (scalar) gluons (D = G) as 

tG c (x,y) = (T X a ' r (x) X b ' s (y)) - (X a ' r (x)} = y) 

iG a {x )V ) = (f X a ' r {x) X h > s {y)) - (X a ' r (x)) (x b ' s (y)) = *G ++ (x,y) 
iG>{x,y) = (x a ' r (x)x b ' s (y)) - (X a ' r (x)) ( X b ' s (y)) =iG + -(x,y) 
iG<(x,y) = (x b ' s (y)x a ' r (x))-(x a > r (x)}(x b ' s (y)) =tG- + (x,y). (A2) 

Here T and T are the usual time and anti-time ordering operators respectively. Our 



convention follows that of Ref. ||12|| . The Wigner Transform of the Green functions 
is given by 

D(X,p) = [d 4 ue tp * u D(x + -,X --) (A3) 



with the center of mass variable X = x + y and the relative distance u = x — y. 

In a standard fashion, the equations of motion for the Wigner transform of these 
Green functions, the so-called transport and constraint equations, can be derived 
0. They read 



- 2ipdxD~ + (X, p) = I coa + I A + I R transport (A4) 



and 



hn x - 2p 2 + 2M 2 ) D~ + (X,p) = I coa + 1^ + 1* constraint, (A5) 

respectively, where M is a generic parton mass, M = m for the gluons and M = 
for the quarks. J co ii is the collision term 

J coll = U- + (X,p)AD + -(X lP ) - n + ~(X lP )AD- + (X,p) 

- rg ain _ floss / A a\ 

— -'coll J coll • K™) 

/-? and l£ are terms containing retarded and advanced components respectively, 
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I R = -U- + (X,p)AD R (X,p) ± D R (X,p)ATl- + (X,p) (A7) 

and 

I A = Il A (X,p)AD- + (X,p) t D- + (X,p)AU A (X,p). (A8) 
In Eqs. (|A6| ) to ( |A8| ) , the operator A is given by 



A := exp (d x ~3 P - t p ~3 x )j (A9) 

and it is set to 1 in the following calculations, which corresponds to the semi-classical 
approximation obtained (h = 1). 

For the calculation of the self-energies, it is useful to introduce the quasiparti- 
cle approximation, in which a free scalar parton of mass M is assigned the Green 
functions 

tD-+(X,p) = I-{6(E p -p )f a (X,p) + 6(E p +p°)f a (X, -p)} (A10) 
iD+-(X,p) = ^{6(E P - p°)f a (X,p) + 6(E p +p°)f a (X, -p)} (All) 

iD-(X,p) = - l„ ; - + 0(-p°)iD+-(X,p) + ®{p Q )iD- + {X,p) 

+ ^r{6(E p - p°)f a (X, p) + 6(E p + p°)f- a (X, -p)} 
&p 

(A12) 



p2 


— M 2 + 


ie 




i 




p'2 


— M 2 + 


ie 




—i 




p2 


— M 2 — 


ie 




—i 




P 2 


— M 2 — 


ie 



iD+ + (X,p) = - — - + e(-p°)iD+-(X,p) + e(p°)iD- + (X,p) 

+ ^{6(E P - p°)f a (X, p) + 6(E P + p°)f- a (X, -p)} 
&p 

(A13) 

with E 2 = p 2 + M 2 , and which are given in terms of the corresponding scalar quark 
and gluon distribution function, f a (X,p), and / a = 1 + f a , where a denotes the 
parton type a = q,g. 
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FIGURES 




FIG. 1. (a) Fock self-energy diagram and (b) Feynman diagram for the process qq — > g 
in lowest order. 




FIG. 2. Two loop self-energy diagrams contributing to £( 2 ) + . 
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